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I review recent results of the high-redshift X-ray selected AGN clustering, based on the 
XMM/2dF survey. Using the luminosity-dependent density evolution luminosity function 
we find that the spatial clustering lengths, derived using Limber's inversion equation, are 
~ 16 and 19 ft -1 Mpc respectively (for the comoving clustering evolution model) while the 
median redshifts of the soft and hard X-ray sources are z ~ 1.2 and 0.75, respectively. Within 
the framework of flat cosmological models we find that these results support a model with 
£l m ~ 0.26, as — 0.75, w~ —0.9 (in excellent agreement with the 3 year WMAP results). We 
also find the present day bias of X-ray AGNs to be b ~ 2. 



Active Galactic Nuclei (AGN) can be detected out to high redshifts and thus their clustering 
properties can provide information on the large scale structure, the underlying matter distribu- 
tion and the evolution with redshift of the AGN phenomenon. From the optical 2QZ and SDSS 
surveys it appears that the QSO clustering properties are comparable to those of local galaxies 
(eg. Croom et al. 2002; 2005; Porciani, Magliocchetti & Norberg 2004; Wake et al. 2004), 
while there is also evidence for the comoving clustering evolution model of active galaxies (see 
also Kundic 1997). Optically selected AGN catalogues however, miss large numbers of dusty 
systems and therefore, provide a biased census of the AGN phenomenon. X-ray surveys, are 
least affected by dust providing an efficient tool for compiling uncensored AGN samples over a 
wide redshift range. From the cosmological point of view an interesting question that remains 
to be addressed is how the high-z X-ray selected AGNs trace the underlying mass distribution 
and whether there are any differences with optically selected samples. 

Early studies of the X-ray AGN clustering, using Einstein and ROSAT data, produced 
contradictory results with other studies finding significant clustering while others not (Boyle & 
Mo 1993; Vikhlinin & Forman 1995; Carrera et al. 1998; Akylas, Georgantopoulos & Plionis 
2000; Mullis et al. 2004). 

Recently, there has been an effort to address this confusing issue and determine the clustering 



properties of both soft and hard X-ray selected AGNs, based on the new XMM and Chandra 
missions (eg. Yang et al. 2003, 2006; Basilakos et al 2004, 2005; Gilli, Daddi, Zamorani 2005; 
Puccetti et al. 2006; Gandhi et al. 2006). Most of these studies find a large correlation length 
for the high-z X-ray AGNs, with hard sources having an even larger correlation length than soft 
sources, with r ~ 17 - 19 hr l Mpc - (eg. Basilakos et al. 2004; Puccetti et al. 2006). 

Here I review our recent results based on the 2 deg 2 XMM/2dF survey, which exploits the 
high sensitivity and the large field-of-view of the XMM-Newton observatory. 

1 The XMM/2df survey: log TV" - log 5 and w(0) 

The XMM-iVeu>fon/2dF survey is a shallow (2-10 ksec per pointing) survey comprising of 18 
XMM-Newton pointings equally split between a Northern and Southern Galactic region near the 
corresponding poles^. Due to elevated particle background we analysed a total of 13 pointings. 
A full description of the data reduction, source detection and flux estimation are presented in 
Georgakakis et al. (2003, 2004). 

We will present results using the soft (0.5-2 keV), hard (2-8 keV) and total (0.5-8 keV) 
band catalogues of the XMM- Newton/ 2dF survey. We only consider sources at off-axis angles 
< 13.5 arcmin. These samples comprise of 432, 171 and 462 sources respectively above the 5a 
detection threshold. The limiting fluxes are /x(0.5 — 2) = 2.7 x 10 -15 erg s _1 cm~ 2 , fx (2 — 8) = 
10~ 14 erg s^ 1 cm~ 2 and /x(0.5 — 8) = 6.0 x 10 -15 erg s _1 cm -2 . 

We derive the source log N — log S after constructing sensitivity maps in order to estimate 
the area of the survey accessible to point sources above a given flux limit (see Basilakos et al. 
2004, 2005). In the 0.5-2 keV band there is good agreement between our results and the Baldi 
et al. (2002) double power-law best fit to the number counts. The Manners et al. (2003) best 
fit is derived for sources in the flux range /(0.5 — 8keV) = 10~ 15 — 8 x 10~ 14 erg s" 1 cm -2 . 
Although our dN/dS is in good agreement with their results in the above flux range, at brighter 
fluxes the surface density of X-ray sources is lower than the extrapolated Manners et al. (2003) 
relation. This suggests that a double power-law is required to fit the 0.5-8 keV dN/dS over the 
flux range io- 15 -io~ 12 erg s 1 cm 2 , which indeed provides an excellent fit to our log N — log S 
(see Basilakos et al. 2005). 

We then calculate the angular correlation function using the estimator: w{9) = /Ndd/Ndr— 
1, of which the uncertainty is: a w = + w(9))/Nor, where N^d and Ndr are the number of 
data-data and data-random pairs, respectively, in the interval [9 — A9, 9 + Ad]. The normaliza- 
tion factor is / = 2Nr/ (Nd — 1), with No and Nr the total number of data and random points, 
respectively. For each XMM pointing we produce 100 Monte Carlo random catalogues having 
the same number of points as the real data which also account for the sensitivity variations 
across the surveyed area (see section 2). Furthermore since the flux threshold for source detec- 
tion depends on the off-axis angle from the center of each of the XMM-Newton pointing, the 
sensitivity maps are used to discard random points in less sensitive areas. This is accomplished 
by assigning a flux to each random point using the source log iV — log S. If that flux is less than 
5 times the local rms noise at the position of the random point (assuming Poisson statistics for 
the background) this is excluded from the random data-set. We have verified that our random 
simulations reproduce both the off-axis sensitivity of the detector as well as the individual field 
log iV — log S. Using the methods described above we estimate w(9) in logarithmic intervals with 
5log9 ~ 0.05. For all three samples we estimate w(9 < 150 ) and find a statistically significant 
signal (see Table 1) at a ^ 3a confidence level (Poisson statistics). We then use a standard x 2 
minimization procedure to fit the measured correlation function assuming a power-law form: 
w(9) = {9 /9) 1 ^ 1 and fixing 7 to 1.8. Note that (a) the fitting is performed for angular sepa- 

a North: RA(J2000) = 13 h 41 m ; Dec.(J2000)=00°00'] and South: RA(J2000)=00 h 57 m , Dec.(J2000)=-28°00' 



X-ray band 


No. of sources 


6* (arcsec) 


X 2 /dof 


P y2 


w(0 < 150") 


0.5-8 keV 


462 


10.8 ± 1.7 


1.50 


0.10 


0.114 ±0.037 


0.5-2 keV 


432 


10.4 ± 1.9 


1.10 


0.35 


0.105 ±0.035 


2-8 keV 


177 


28 ±9 


0.88 


0.57 


0.128±0.080 



Table 1: Angular correlation function analysis results with their la (A\ 2 = 1.00) uncertainties. The fits are 

produced after imposing 7 = 1.8. 
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Figure 1: The XMM/2dF survey two-point angular correlation function for the 3 bands considered. 



rations in the range 40-1000 arcsec, and (b) our results are insensitive to both the upper cutoff 
limit in 9 and the angular binning (for more than 10 bins) used to estimate w(9). The resulting 
raw values of o are corrected for the integral constraint and the amplification bias (see Vikhlinin 
& Forman 1995 and Basilakos et al. 2005), although such corrections are quite small. The final 
results are presented in Table 1. 



2 The spatial correlation length of the XMM/2dF soft X-ray sources 

The spatial correlation function can be modeled as (eg. de Zotti et al. 1990): 

£(r,z) = (r/r )-~< x (I + z)~^ , (1) 

where e parametrizes the type of clustering evolution. For e = 7 — 3 (ie., e = —1.2 for 7 = 1.8), 
the clustering is constant in comoving coordinates (comoving clustering), a model which appears 
to be appropriate for active galaxies (eg. Kundic 1997; Croom et al. 2005). 

In order to invert the angular correlation function to three dimensions we utilize Limber's 
integral equation (eg. Peebles 1993). For a spatially flat Universe, Limber equation can be 
written as: 

Jo° Io° x 4 4> 2 {x)£{r, z)dxdu 
[f °° x^(x)dx]i 

where 4>(x) is the distance selection function (the probability that a source at a distance x is 
detected in the survey), x is the proper distance related to the redshift through (see Peebles 
1993): 

x{z) = — jf — - with E{z) = [Ml + z? + ^a] 1/2 , (3) 
The selection function, 4>(x), is related to the number of objects in the given survey with a solid 
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Figure 2: Left Panel: The expected redshift distribution (continuous line) by the Miyaji et al (2000) LDDE 
soft-band luminosity function. Right Panel: The corresponding redshift distribution for the Ueda et al. (2003) 
LDDE hard-band luminosity function (the dashed line is based on the La Franca et al. 2005 luminosity function) . 
The histogram in the left panel corresponds to the distribution of the Schmidt et al (1998) X-ray sources of the 
ROSAT Lochman Deep Field (albeit having a flux limit slightly lower than of our survey), while that of the right 
panel on limited spectroscopic and photo-z data of our XMM/2dF survey. 



angle O s and within the shell (z, z + dz), by: 

^L = n sX ^( x )^E- 1 (z) . (4) 

Since we do not have complete redshift information for our sources we estimate dN/dz using 
the X-ray source luminosity function and the specific flux-limit of our samples, via the relation: 
4>{x) = . Q(L X , z)dL, where $>(L x ,z) is the luminosity dependent density evolution lumi- 
nosity (LDDE) function. For the soft-band we use that of Miyaji, Hasinger & Schmidt (2000) 
while for the hard-band that of Ueda et al (2003). In Fig. 2 we present the expected redshift 
distributions of the soft and hard X-ray sources together with the histogram of some limited 
spectroscopic and photo-z data (see caption for details). The LDDE model predicts a redshift 
distribution with a median redshift of z ~ 1.2 and 0.75 for the soft and hard sources respectively. 
Finally, the expression for w(6) satisfies the form: 




Note that, the physical separation between two sources, separated by an angle 9 considering the 

1 /2 

small angle approximation, is given by: r ~ (1 + z)^ 1 (v? + x 2 9 2 ) 

Then the inversion of eq. (5), using the LDDE luminosity evolution model, e = —1.2 and 
the concordance cosmological model, provides a spatial correlation length of r ~ 16.4 ± 1.3 h~ l 
Mpc and ~ 19 ± 1.3 h^ 1 Mpc, for the soft and hard bands, respectively. These results are in 
very good agreement with a recent XMM based study of the ELAIS-S1 field by Puccetti et al. 
(2006) and comparable to those of Extremely Red Objects (EROs), of luminous radio sources 
(Roche, Dunlop & Almaini 2003; Overzier et al. 2003; Rottgering et al. 2003) and of bright 
distant red galaxies (Foucaud et al. 2006). 

Our r values, however, are significantly larger than those derived from optical AGN surveys 
(which trace mostly the unobscured component): r a ~ 5.4 — 8.6 h^ 1 Mpc (eg. Croom & Shanks 
1996; La Franca et al. 1998; Croom et al. 2002; Grazian et al. 2004; Porciani et al. 2004; 
Wake et al. 2004). We can push our inverted r values to approximate closely the optical AGN 



results only if we use the constant in physical coordinates clustering evolution model (e = —3), 
in which case we obtain r Q ~ 7.5 ± 0.6 Yr x and ~ 13.5 ± 3 hr x Mpc (for the soft and hard bands 
respectively). 

3 Cosmological Constraints 

It is well known (Kaiser 1984) that according to linear biasing the correlation function of the 
mass-tracer (£ bj) and dark- matter one (£dm)> are related by: 

£ob](r,z) =b 2 (z)£ DM (r,z) , (6) 

where b(z) is the bias evolution function. In this study we use the bias model of Basilakos & 
Plionis 2001; 2003) which is based on linear perturbation theory and the Friedmann-Lemaitre 
solutions of the cosmological field equations. We quantify the underlying matter distribution 
clustering by presenting the spatial correlation function of the mass ^dm^z) as the Fourier 
transform of the spatial power spectrum P(k): 

&M(r ,^(i±|i^/;, 2m ^«, (7) 

where k is the comoving wavenumber and e = —1.2, according to the constant in comoving 
coordinates clustering evolution model. As for the power spectrum, we consider that of CDM 
models, where P{k) = Pok n T 2 (k) with scale-invariant (n = 1) primeval inflationary fluctuations. 
In particular, we use the transfer function parameterization as in Bardeen et al. (1986), with the 
corrections given approximately by Sugiyama (1995). The normalization of the power spectrum 
is given by: 

(8) 

where as is the rms mass fluctuation on R = 8/i _1 Mpc scales and W(kR) is the window function. 
Note that we also use the non- linear corrections introduced by Peacock & Dodds (1994). 

We have chosen to use either the standard normalization given by: as — 0.50"^ with 
7 ~ 0.21 — 0.22w + 0.33Sl m (Wang & Steinhardt 1998), or to leave as a free parameter to be 
fitted by our analysis. 

3.1 X-ray AGN Clustering likelihood 

It has been shown that the application of the correlation function analysis on samples of high 
redshift galaxies can be used as a useful tool for cosmological studies (eg. Matsubara 2004). In 
what follows we review a similar analysis, presented in Basilakos & Plionis (2005, 2006), utilizing 
a x 2 likelihood procedure to compare the measured XMM soft source angular correlation function 
with the prediction of different spatially flat cosmological models. In particular, we define the 
likelihood estimator as: £ AGN (c) oc exp[— Xagn( c )/2] with: 

n 

Xagn( c ) = J2 
1=1 

where c is a vector containing the cosmological parameters that we want to fit and ai the 
observed angular correlation function uncertainty. We assume a flat (r2 tot = 1) cosmology with 
primordial adiabatic fluctuations and baryonic density of ^l-^h? ~ 0.022 (eg. Kirkman et al. 
2003; Spergel et al. 2006). In this case the corresponding vector is c = (fJ m , w, as, h, b Q ) and we 
sample the various parameters as follows: the matter density fi m € [0.01, 1] in steps of 0.01; the 
equation of state parameter w G [—3, —0.35] in steps of 0.05, the power spectrum normalization 
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Figure 3: Likelihood contours in the (w, ft) plane (left panel), the (o"8,f2 m ) plane (central panel) and the (u>,f2 m ) 
plane (right panel). The contours are plotted where — 21n£/£ max is equal to 2.30, 6.16 and 11.83, corresponding 
to la, 2a and 3a confidence level. In the (S7 m , w) plane we plot as thick lines the likelihood contours derived from 

the SNIa Hubble relation. 



o"8 € [0.4,1.4] in steps of 0.02, the dimensionless Hubble constant h S [0.5,0.9] in steps of 0.02 
and the X-ray sources bias at the present time b E [0.5,4] in steps of 0.05. Note that in order 
to investigate possible equations of state, we have allowed the parameter w to take values below 
-1. Such models correspond to the so called phantom cosmologies (eg. Caldwell 2002). 

The resulting best fit parameters for the e = —1.2 clustering evolution model are presented 
in Table 2. In the first two rows we present results based on the traditional Wang & Steinhardt 
(1998) as normalization. Note that our estimate of the Hubble parameter h (left panel in Fig. 3) 
is in very good agreement with those derived (h = 0.72 ±0.07) by the HST key project (Freeman 
et al. 2001). In the last two rows of Table 2 we leave as free but fix the Hubble constant to 
h = 0.72. In Fig. 3 we present the la, 2a and 3a confidence levels in the various parameter 
planes, marginalizing over the rest of the parameters. We find that w is degenerate with respect 
to both h and the bias at the present time. 

When we leave free the as parameter, our fit (central panel of Fig. 3) provides a value which 
is in excellent agreement with that derived by the recent 3-years WMAP results (Spergel et al. 
2006) b . Therefore, allowing for the first time values w < — 1 (Phantom models) we can derive 
a (r2 m ,o"8) relation, a good fit of which is provided by : 

a 8 = 0.34(±0.01) n m ^ Q ^ (10) 

with j(n m ,w) = 0.22(±0.04) - 0.40(±0.05)w - 0.052(±0.040)S! m . 

Note that eq. (10) produces as values which are significantly smaller than the usual cluster 
normalization (Wang & Steinhardt 1998) but are in good agreement with the 3-years WMAP 
results; for example for w ~ — 1 and Q m ~ 0.28 we get as ^ 0.73 ± 0.03. 

Inspecting the thin contours in the right panel of Fig. 3 it becomes evident that w is 
degenerate, within the la uncertainty, with respect to Q m . Therefore, in order to put further 
constraints on w we additionally use a sample of 172 supernovae SNIa (see Tonry et al. 2003). 

3.2 The AGN+SNIa likelihoods 

We combine the X-ray AGN clustering properties with the SNIa data by performing a joined 
likelihood analysis and marginalizing the X-ray clustering results over as, h and 6q- The vector 
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Table 2: Cosmological parameters from the likelihood analysis. Errors of the fitted parameters represent la 
uncertainties. Note that for the joined analysis (2nd and 4th rows) the corresponding results are marginalized 
over the parameters that do not have errorbars, for which we use the values indicated. 



c now becomes: c = (Q m , w). The SNIa likelihood function can be written as: £ SNIa (c) oc 

ex P[-XSNIa( c )/ 2 ]> With: 
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where Dl(z) is the dimensionless luminosity distance, Dl(z) = H (l + z)x(z) and z% is the 
observed redshift. The results are shown as thick lines in Fig. 3 and represent the la, 2a, and 
3a, confidence levels. The joint likelihood function (D omt (Q m , w) = £ AGN x £ SNIa ) peaks at: 
Sl m = 0.26 ± 0.04 with w = — 0.90^o'o 5 . Using eq. (10) we find that the normalization of the 
power spectrum that corresponds to these cosmological parameters is ag ~ 0.73. It should be 
pointed out that our results are in excellent agreement with the recent 3-year WMAP results of 
Spergel et al. (2006). 

Other recent analyzes, utilizing different combinations of data, seem to agree with our results. 
For example, Sanchez et al. (2006) used the WMAP (1-year) CMB anisotropies in combina- 
tion with the 2dFGRS power spectrum and found £l m ~ 0.24 and w ~ —0.85, while Wang & 
Mukherjee (2006) utilizing the 3-years WMAP data together with SNIa and galaxy clustering 
results found w ~ —0.9. 
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